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The role of monopoles in the deconfinement transition is discussed in the 
framework of abelian projection in the maximally abelian gauge in T ^ 
SU{2) QCD. Only one (or a few near (3 C ) long connected monopole loop 
exists uniformly through the whole lattice in each vacuum configuration in 
addition to some very short loops in the confinement phase and the long 
loop disappears in the deep deconfinement region. Energy-entropy balance 
of the long loops of i 
of the deconfinemen 
critical temperature. 
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of the long loops of maximally extended monopoles explains the existence 
of the deconfinement transition and reproduces roughly the value of the 



§1. Introduction 

Modern computers enable us to simulate the system of quarks and gluons. It has 
been evident from Monte-Carlo simulations of lattice QCD that quarks are confined. 
However why and how quarks are confined is not yet known. To understand the mech- 
anism of the confinement is important in order to explain hadron physics out of QCD. 
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The 'tHooft idea of abelian projection of QCD is interesting.!^ The abelian projection 
is to fix the gauge in such a way that the maximal torus group remains unbroken. Af- 
ter the abelian projection, monopoles appear as a topological quantity in the residual 
abelian channel. QCD can be regarded as an abelian theory with electric charges and 
monopoles. If the monopoles make Bose condensation, charged quarks and gluons are 
confined due to the dual Meissner effect. 

There are, however, infinite ways of extracting such an abelian theory out of 
SU(3) QCD. It seems important to find a good gauge in order to test whether the 
conjecture is true or not at least on a practically accessible lattice. It has been found 
that, if one adopts a gauge called a maximally abelian (MA) gauge^'cPo the 't Hooft 
conjecture is seen to be beautifully realized in QCD. 

1. The first interesting findings are some phenomena called abelian dominance. 
Abelian loop operators composed of abelian link variables alone seem to re- 
produce essential features of color confinement in the MA gauge .EP® Explicitly, 
abelian static potentials derived from abelian Wilson loops can reproduce the 
string tensions which is a key quantity of confinement. Abelian Polyakov loops 
and abelian energy densities play the role of an order parameter of the deconfine- 
ment transition in finite-temperature pure QCD.lf The abelian quantities show 
clearer behaviors around the critical coupling. 

2. Furthermore, it has been shown that monopoles alone are responsible for the 
above interesting phenomena played by the abelian quantities. Abelian Wilson 
loops^'Q^tP'tlJ) and abelian Polyakov loopsEJ^iiHP are written by a product of 
monopole and photon contributions. The string tension and the behavior of the 
Polyakov loop as an order parameter of the deconfinement transition are all re- 
produced by the monopole contributions alone. 

3. Abelian dominance suggests that there must exist a U(l) invariant effective ac- 
tion which can explain confinement. It is possible to get the U(l) effective action 
in terms of monopole currents in T = SU(2) QCD after performing a dual 
transformation numerically.^^ a) c§>0 The effective action is fixed also for ex- 
tended monopolesHP Considering extended monopoles corresponds to making a 
block spin transformation on the dual lattice. The actions determined are very 
interesting because they appear to satisfy a scaling behavior of the renormalized 
trajectory on which one can take the continuum limit Ji3 IZPO) \}3 If the similar be- 
haviors would show up on a larger lattice, one could conclude that T = SU{2) 
QCD is always ( for all j3 in the infinite lattice limit) in the monopole condensed 
phase. Color confinement would be proved then. 

What happens in the case of the deconfinement transition in the finite-temperature 
(T 7^ 0) QCD? Can the similar monopole dynamics like the energy-entropy balance 
explain the transition? This is the main subject of this note. It is highly probable, since 



the monopole contributions alone explain the string tension and the behavior of the 
Polyakov loop also in T ^ QCD.B©0© 

In the next section, what is abelian projection is shortly reviewd. Section 3 reviews 
the results in T = QCD. Section 4 is the most important part of this note discussing 
the case of T ^ SU(2) QCD. Final section is devoted to conclusion and remarks. 

§2. Abelian projection 

2.1. Abelian projection in the continuum QCD 

Abelian projection of QCD is done as follows. Choose an operator X(x) which is 
transformed non-trivially under SU(3) transformation: 

A^x) -» A fl (x) = V(x)A fl (x)V^(x)--d fl V(x)V\x) (1) 

ip(x) -► ${x) = V(x)ip(x). (2) 

Abelian projection is to choose V(x) so that X(x) is diagonalized as 

( Xx{x) \ 

X(x)^X(x)=\ X 2 (x) " . (3) 

V Ux)J 

V(x) is fixed up to the diagonal element d(x) of SU(3), where 

/ exp(iai(a;)) Q \ 3 

d(x)=\ exp(z« 2 (x)) ~ 6^(3), ^a,(x) = 0. (4) 

\ U exp(ia 3 (x)) / i=1 

{d(x)} is the maximum torus group of SU(3), which is the residual U(l) x U(l) gauge 
symmetry. 

We look at QCD at this stage without further fixing the gauge of the residual 
symmetry. First, we explore how the fields after the abelian projection transform under 
an arbitrary SU(3) gauge transformation W(x). Since X(x) is a functional of (gauge) 
fields and so it transforms under W(x). Hence V(x) also transforms non-trivially. Let 
us fix the form of V(x) such that all diagonal components of the exponent of V(x) are 
zero. This is always possible if one uses the residual symmetry. Then V(x) is found to 
transform under W(x) as 

V{x) -^ V w (x) = d w (x)V{x)W\x). (5) 

V w (x) diagonalizes an operator which is transformed from X(x) under W(x). d w (x) 
is necessary for V (x) to take the fixed form. 



The gauge field after the abelian projection, A p (x), transforms under W(x) as 
A^x) ^U A^(x) = d w {x)A^x)d w \x) - -d^d w {x)d w \x). (6) 
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After the abelian projection, A p (x) transforms only under the diagonal matrix d w (x). 
Since the last term of (||) is composed of the diagonal part alone, the diagonal part 
of A^x) transforms like a photon. The off diagonal part of A p (x) transforms like a 
charged matter. The quark field transforms under W(x) as 

4>{x) -^ $ w (x) = d w (x)ij(x). (7) 

It is important that ipj (x)ipi(x) and ^i^^ix)^)^) are neutral and at the same time 
invariant under any SU(3) transformation W(x). Color confinement is regarded as 
abelian charge confinement after abelian projection. 

The most interesting fact of abelian projection is that monopoles appear in the 
residual abelian channel. We treat SU(2) QCD for simplicity. After abelian projection, 
we define an abelian field strength as 

f^(x) = dffix)-dffix). (8) 

The abelian field A 3 (x) written in terms of the original field is 

Al(x) = Y%x)A;(x) - - I e 3ab Y a (x)d,Y b (x), (9) 

9 1 + Y 6 {x) 

where Y(x) = V\x)a^V (x) = Y a (x)o a . Y a (x) obeys 

Y a {x)Y a (x) = 1. (10) 

fnu{x) can be rewritten in the form 

UAx) = d^Y a (x)A a v (x)) - d v (Y a {x)Afa)) - -e abc Y a {x)d l j\x)d u Y c {x) (11) 

in terms of the original field. A current 

K(x) = \e, vpa d v r{x) (12) 

= ^e liVpa e abc d"Y'\x)d p Y\x)d"Y c {x) (13) 

^9 

is always zero if V(x) is fixed. However, at a point x where the eigenvalue of the 
diagonalized operator X(x) is degenerate, V(x) is not well defined and k p {x) does not 



vanish there. We calculate the charge in the three dimensional volume f2 around xxB 

g m = I h(x)d 3 x = ±- f e 0upa e abc d»Y a (x)d p Y\x)d°Y c (x)d 3 x (14) 

Jn 2g Jn 

= 7T [ e ijk e abc Y a (x)d j Y b (x)d k Y c (x)d 2 a i (15) 

Zg JdQ 
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(16) 



where n is an integer, n is a topological number corresponding to a mapping between 
the sphere (0) in the parameter space and the sphere dQ of Q. Because this equation 
represents the Dirac quantization condition, g m can be interpreted as a magnetic charge. 
The monopole current fc M (x) is a topologically conserved current d fl k fl (x) = 0. Abelian 
projected QCD can be regarded as an abelian theory with electric charges and monopoles. 
'tHoofuiP conjectured if the monopoles condense, abelian charges are confined due to 
the dual Meissner effect. This means color confinement. 

2.2. Abelian projection on a lattice 

We can perform abelian projection on a lattice similarly. Choose an operator 
X(s) in, for simplicity, 577(2) QCD. The gauge transformation on a lattice is 

U(s,fi)=V(s)U(s,fi)V\s + fi), (17) 

where U(s, ft) represents a link field corresponding to a gauge field in the continuum 
theory. After abelian projection is over, abelian link fields can be separated from 577(2) 
link fields as follows: 



= , y i-|«W -^>)^f o (18) 



c(s,n) V 1 — |c(s,/x)| 



= C(s,p)u(s,p). (19) 

u(s,fi) transforms like a photon and C(s,fi) transforms like a charged matter under 
the residual U(l) gauge symmetry. The abelian field strength is defined as a plaquette 
variable 

W s ) = ^( s ) + e As + A) - O^s + i>) - 9 v {s). (20) 

The monopole current is defined^ as 

k v( s ) = ifnvpad v n p(J (s + ft), (21) 

where d v is a forward derivative on a lattice and 9 pu (s) is decomposed into 

M s ) = M s ) + 2t™^(s), -tt < 6 pu (s) < tx. (22) 



n^ u (s) is an integer corresponding to the number of the Dirac string through the pla- 
quette. The monopole currents are conserved topologically 

d»K = °» ( 23 ) 

where d' is a backward derivative on a lattice. The monopole currents make closed 
loops on the four dimensional lattice. The above monopole current is defined by n^ v {s) 




Fig. 1. An extended cube on which an 2 3 extended monopole is defined as the sum of eight 
(=2 3 ) smallest monopoles. 



surrounding the smallest cube as shown in fl21|) . To study the lpne range behavior 
important in QCD, considering extended monopoles^ is essential .H^ollM'llB) They are 
defined by n At!/ (s) surrounding an extended cube. For example, Fig.l represents a 2 3 
cube defining a 2 3 extended monopole. Adopting a n 3 extended monopole corresponds 
to performing a block spin transformation on a dual latticeo^O'il^ and so is suitable 
for exploring the long range property of QCD. When adopting a n 3 extended monopole 
on A^ 3 x Nt lattice, the effective lattice on which the extended monopole runs is 



A", 



n 



N t 



n 



(24) 



which we call a renormalized lattice. 



§3. Monopole condensation in T = QCD 

As shown in the introduction, a gauge called a maximally abelian (MA) gauge 
is very interesting. In the MA gauge, 

x ( s ) = H [ u ( s : A)o"3^ f (s, A) + U j (s - fa, p,)a 3 U(s - jx, p) 



(25) 
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Fig. 2. Coupling constants /i versus j3 for 2 3 ,3 3 , and 4 3 extended monopoles on 24 4 lattice. 

is diagonalized. Abelian loop operators composed of u(s, fi) alone seem to reproduce 
essential features of color confinement. Here we review briefly the results showing that 
monopoles are a key quantity of confinement in T = SU{T) QCD. They are obtained 
by Shiba and one of the authors (TS) recentlySi&t'Boi]) 

The first interesting result is that one can determine an effective U(l) action in 
terms of monopole currents in the MA gauge in SU(2) QCDJI^'lPH 4 !'^ The partition 
function of interacting monopole currents is expressed as 



Z = (U E )(nW)>°) ex p(-^D- 



(26) 



s iA* k li (s)=— oo 



It is natural to assume S[k] = J2i fiSi[k]. Here fi is a coupling constant of an in- 
teraction Si[k]. For example, /i is the coupling of the self energy term En l( j(^( 5 )) 2 ) 
fi is the coupling of a nearest-neighbor interaction term J2 n ,ti k fl (s)k IJi (s + «land f% 
is the coupling of another nearest-neighbor term Ylm.^v kfi(s)kJs + t> ) .113* 'E^P Shiba 
and Suzukill^oilflO extended a method developed by SwendsenEP to the system of 
monopole currents obeying the current conservation rule (|23|) . The monopole actions 
are obtained locally enough for all extended monopoles considered even in the scaling 
region. They are lattice volume independent. The coupling constant f\ of the self- 
energy term is dominant and the coupling constants decrease rapidly as the distance 
between the two monopole currents increases. 

Since the action is fixed, it is possible to study energy and entropy balance of 
monopole loops in order to confirm the occurence of monopole condensation. If the 
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Fig. 3. Coupling constants f\ versus b. 



entropy of a monopole loop exceeds the energy, the condensation of a monopole loop 
occurs. As done in compact QED,clP the entropy of a monopole loop can be estimated 
as In 7 per unit loop length. Since monopole currents are distributed randomly in 
average for large L, interaction terms between two separate currents is expected to be 
canceled. That this actually happens will be shown later in the case of T ^ QCD case. 
Hence the action may be approximated by the self energy part fiL. Here dominance 
of currents with a unit charge is used. Since /i is regarded as the self energy per unit 
monopole loop length, the free energy per unit monopole loop length is approximated 
by 



A -In 7. 



(27) 



If /i < In 7, the entropy dominates over the energy, which means condensation of 
monopoles. In Fig. 2, f\ versus (3 for various extended monopoles on 24 4 lattice is 
shown in comparison with the entropy value In 7. Each extended monopole has its own 
j3 region where the condition fi < In 7 is satisfied. When the extendedness is bigger, 
larger /3 is included in such a region. Larger extended monopoles are more important 
in determining the phase transition point. 

The behaviors of the coupling constants are different for different extended mono- 
poles. However, if we plot them versus b = n x a(/3), we get a unique curve as in Fig.||. 
The coupling constants seem to depend only on b, not on the extendedness nor f3. There 
is a critical b c corresponding to critical /?", i.e., b c = na(j3™). The monopole action may 



be fitted by 

— 1 Att 

S[k] = T, m obk,(s)k,(s) + -(^) 2 Y,k,(s)D(s-s')k,(s% (28) 

where g(b) is the SU(2) running coupling constant 

^"=2^ ln (^) + 4i^ hUn (^- (29 > 

D(s) is a modified lattice Coulomb propagator. This form of the action is predicted 
theoretically by Smit and Sijs.c^P The solid line is the prediction given by the action 
with the parameters written in Fig.^. 

Suppose the effective monopole action remains the same for any extended mono- 
poles larger than 4 3 in the infinite volume limit. Then the finiteness of b c = na(/3™) 
suggests (3™ becomes infinite when the extendedness n goes to infinity. SU(2) lattice 
QCD is always (for all /3) in the monopole condensed and then in the color confinement 
phased This is one of what one wants to prove in the framework of lattice QCD. 

Notice again that considering extended monopoles correponds to performing a 
block spin transformation on the dual lattice. The above fact that the effective actions 
for all extended monopoles considered are the same for fixed b means that the action 
may be the renormalized trajectory on which one can take the continuum limit. Our 
results suggest the continuum monopole action takes the form (p8|) predicted by Smit 
and Sijs.E^ The simulation of the monopole action is in progress. 

Shiba and SuzukffCPll) showed furthermore that monopoles alone can reproduce 
the full value of the string tension in SU(2) QCD. These facts strongly suggest that 
the monopole condensation occurs in QCD and quarks and gluons are confined due to 
the dual Meissner effect. 

§4. Deconfinement transition in T^O QCD 

Now let us study the role of the monopoles in the deconfinement transition in 
T^O QCD. There have been interesting data already suggesting the importance of 
monopoles also inT ^ QCD. Monopole contribution to the string tension has been 
studied in the MA gauge also in T ^ SU(2) QCD. The string tension from the 
monopoles almost agrees with that from the abelian Wilson loops in the confinement 
phase, whereas it vanishes clearly in the deconfinement phase. The data are shown in 
Fig.4. The string tension from the photon part is negligibly small. 

An abelian Polyakov loop which is written in terms of abelian link fields alone is 
given by a product of contributions from Dirac strings of monopoles and from photons. 
It has been foundtiPO that the characteristic features of the Polyakov loops as an order 
parameter of the deconfinement transition are due to the Dirac string contributions 
alone. In Fig. 5, the SU(2) data in the MA gauge are plotted. The abelian Polyakov 
loops vanish in the confinement phase whereas they begin to rise for (5 larger than the 
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Fig. 4. Physical string tensions (circle) and spatial string tensions (square) from monopolcs 
on 24 3 x 8 lattice. 
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Fig. 5. Monopole Dirac string and photon contributions to Polyakov loops in the MA gauge 
in SU(2) QCD on 16 3 x 4 lattice. Total means the abelian Polyakov loops. 
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critical temperature f3 c = 2.298. The Dirac string contribution shows similar behaviors 
more drastically. The photon part has a finite contribution for both phases and it 
changes only slightly. The fact that monopoles are responsible for the essential feature 
of the Polyakov loop is found also in U(l) and in SU(3) in the MA gauge. In addition, 
a remarkable result has been found. The behavior of the abelian Polyakov loops as 
an order parameter and the monopole responsibility are seen also in other gauges like 
unitary ones. This is the first phenomena suggesting gauge independence of the 'tHooft 
conjecture. 

4-1. The expected role of maximally extended monopoles 

In T = QCD, each extended monopole has its own j3 region where the entropy 
dominance of monopoles occurs, i.e., where the condition f\ < In 7 is satisfied. When the 
extendedness is bigger, larger /3 is included in such a region. Larger extended monopoles 
are more important in determining the phase transition point. In the infinite- volume 
limit, one can adopt infinitely extended monopoles. If the situations remain unaltered 
even in the infinite- volume limit, one can prove that T = SU(2) lattice QCD is always 
(for all j3) in the monopole condensed and then in the color confinement phase.a-* 

What happens inT ^ QCD? It is expected that the confinement - deconfine- 
ment transition in T ^ QCD also can be explained in terms of monopoles. Then 
energy and entropy balance of monopoles must be the mechanism as shown above in 
T = QCD. However, there is a big difference between T = and T ^ QCD. 
Namely the time extent N t is kept finite in T 7^ case. This means that the extended- 
ness of monopoles running in the space direction (which we call dynamical monopoles) 
is restricted to be finite, since monopoles are defined by a three dimensional cube per- 
pendicular to the direction of the current. The monopole currents which contribute 
to the physical string tension are those running in the direction perpendicular to the 
Wilson loops. Because the Wilson loops are on a plane which include the time axis and 
a space one, dynamical monopoles are essential in the confinement mechanism. It is our 
expectation that the existence of the maximum size of extended dynamical monopoles 
is the reason there is a deconfinement phase transition at a finite value of (3 C , i.e., T c 
in the finite temperature QCD. Here T c is the critical temperature. 

In the following subsections, we explore what is the maximum size of extended 
monopoles and whether /3 C can be explained by balancing of energy and entropy of the 
maximally extended monopole loops in the MA gauge. 

4-2. Monopole action in T 7^ QCD 

Monopole action can be calculated from monopole current configurations sim- 
ilarly as in T = SU(2) QCD. The partition function is expressed by Eg. fl2"S|) . In 
contrast with the case of T = QCD, we must treat time and space directions sep- 
arately. Since the number of sites in the time direction is small, we consider only up 
to next to the nearest quadratic couplings. Owing to the current conservation, five 
interaction terms are sufficient in this case. In Fig. 6 the action is calculated in both the 
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Fig. 6. Coupling constants /i(circle), /^(square) and ^(diamond) versus /3 on 16 3 x 4 lattice 
for l 3 monopoles. Open symbols are the couplings connecting two space-like currents and 
filled ones are those connecting two time-like currents. 

confinement and the deconfinement phases on 16 3 x 4 lattices. We used 50 (100 at some 
points) configurations adopting l 3 monopoles. /i is the coupling of the self energy term 
En,/i(^( 5 )) 2 ) fi is the coupling of the term J2 n ,u k ll (s)k tl (s + (x) and f 3 is the coupling 
of the term, En,u^^(*)^( s + ^)- Other couplings between currents larger distance 
apart are much smaller, fi is dominant in both phases. In the deconfinement phase, 
the discrepancy between space-space and time-time couplings is large, whereas it is 
negligible in the confinement phase. 

/i of the monopole action is plotted in Fig. 7 in the confinement phase on various 
lattices. In the confinement phase the action is independent of N s and N t ,i.e., of the 
lattice volume and the action obtained seems the same as that given in T = QCD. 

Now that the energy of a monopole loop is evaluated, let us check if all terms 
other than the self energy are canceled because of randomness of the monopole loop 
distribution. In Fig. 8 we compare the value of the total monopole action with that of 
the self energy term and also with that of f± x L, where L is the length of a monopole 
loop. The length of monopole loops will be defined in detail in the next subsection. 
The data show that the cancellation is good and the total action is well approximated 
only by the self energy term. Furthermore the value of the self energy term is in fairly 
good agreement with fi x L, which suggests dominance of monopole currents with unit 
charge. Next, we pay attention to the entropy of a monopole loop. 
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Fig. 7. Coupling constants f\ in the confinement phase versus /3 on various lattices. The 
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Table 1. The total number of loops of length L on a sample of 10 configurations on 16 3 x 4 
lattice at (5 =2.20 for l 3 monopoles (a) and 2 3 monopoles (b). 

4-3. Monopole loop length 

Suppose that the behavior of a monopole loop is represented by non-backtracking 
random walks, the number of loops of length L on four- dimensional large lattice behaves 
as 7 L . The entropy of a monopole loop can be estimated as In 7 per unit loop length 
in T = QCD. However, the entropy of the loop inT ^ QCD must be affected by 
the finiteness of the time direction and the entropy may become smaller than In 7. To 
estimate the entropy, we have to investigate the behavior of monopole loops carefully 
in T ^ QCD. 

First, we calculate the length of monopole loops. Bode et al.EP studied the length 
of monopole loops in U(l) lattice gauge theory. They found that in the confinement 
phase one large loop occupies the whole lattice and in the Coulomb phase it splits into 
several smaller pieces. We define the monopole loop length following Bode et && in 
the MA gauge in SU(2) QCD. If two closed loops occupy some sites in common, these 
two loops are regarded as one loop with the number of crossings. It is found that in 
the confinement phase, a long monopole loop exists in each configuration in addition 
to some short loops. The separation of long loop and other short loops is clearly seen 
in Table 1, where the data are obtained from 10 configurations at /3 = 2.20 on 16 3 x 4 
lattice. Note that the number of long loop is equal to the number of configurations. In 



14 



this case, f3 c is 2.298. For (3 < j3 c long monopole loops have a characteristic length at 
each j3 and the long loops become shorter or split into a few parts around f3 c . In the 
deep deconfmement region, no long loop exists and all monopole loops are short. These 
results are consistent with those in the U(l) case.c^ 
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Table 2. The length of the longest loop, the number of crossings and the ratio (percent) of the 
number of multicharged monopole currents to the total number of monopole currents in the 
case of 2 3 -extended monopoles (a) and that of 4 3 -extended monopoles (b) on 24 3 x 8 lattice. 
The numbers are averages over 30 configurations. 

In Table 2, the /3-dependence of the loop length (L), the number of crossings and 
the ratio of the number of multicharged currents to the number of the total currents are 
shown. Those quantities are calculated from 30 configurations for two different sized 
monopoles on 24 3 x 8 lattice. /3 C is 2.51 on this lattice. The number of crossings is 
small in the case of 2 3 monopole, whereas it is not so small in the case of 4 3 monopole. 
The monopole currents with a unit charge are dominant for both extended monopoles, 
which is consistent with that f\ is dominant on the monopole action. In the case of 4 3 
monopole, (L) exceeds the site number of the renormarized lattice for j3 < /3 C , whereas 
they are nearly equal around /3 C . In the case of 2 3 monopole (L) is smaller than the 
number of the sites. 

Recently, Ejiri et alJa'a) have evaluated the string tension derived from long 
and short monopole loops separately. They have found that the long loops alone can 
reproduce the full value of the string tension in spite of the fact that only a few 
percents of a total dual links are occupied by the long loop. Other short loops do not 
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Fig. 9. R 2 versus (5 on various lattices. The number in the parentheses represents the extend- 
edness of monopoles. 



contribute to the string tension. This suggests that the long loops are essential in the 
deconfinement transition. 

4-4- The distribution of monopole currents 

Now we pay attention to long loops and investigate the distribution of monopole 
currents contained in a long loop. Even the longest loop occupies only about 10 percent 
of the total dual links. Does the long loop spread uniformly through the lattice or not? 



In Ref.^3|), it is filling up the whole lattice in compact QED. 

We calculate the mean square of the distance from the center as follows: 



R 2 



1 



£( 



Raf 



(30) 



where L is the length of the monopole loop, r*j is the position of a monopole current 
and Rq is the position of the center. If the monopole currents exist uniformly on the 
N 3 x N t lattice, R 2 should be evaluated as 



N 2 N 2 



f N a /2 f N 3 /2 f N 3 /2 f N t /2 jy 2 

/ dx dy dz dt(x 2 + y 2 + z 2 + t 2 ) = — + 

J-N s /2 J-N./2 J-N./2 J-Nt/2 4 



NfN t J-N s /2 J-N s /2 ' J-N s /2 J-Nt/2 

In the case of n 3 extended monopole, R 2 should be 

(N s /n) 2 + (N t /n) 2 



12 



1 



12 



(31) 



(32) 
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Our data of R 2 are shown in Fig.9. f3 c is 2.426 in the case of N t = 6 and (3 C is 2.51 
in the case of Nt — 8. In the confinement phase, the data almost coincide with the 
expected values (p2|) . For example, R 2 is expected to be 147 in the case of l 3 monopole 
on 24 3 x 6 lattice and in the case of 3 3 monopole R 2 is expected to be 16. The long 
monopole loop is almost uniformly distributed in the whole lattice in the confinement 
phase. In the deconfinement phase R 2 starts to decrease rapidly above f3 c . In the case 
of extended monopoles the uniformity is better than in the l 3 case for (3 < (3 C . 

4-5. Monopole effective size 

Assuming that crossings of the loop are ignored, an n 3 monopole loop seems to 
spread almost uniformly through the whole lattice for f3 < (3 C . Considering that almost 
all magnetic charge is ±1 and that the loop length is rather restricted in comparison 
with all link numbers 4 x N 3 x N t , there appears to work some repulsive force between 
currents. To study the force, we may define an effective size of monopoles as follows: 

,(„) = W*W/"> , (33) 

where the numerator is the volume of the renormalized lattice of n 3 extended monopoles 
and the denominator is the length of the monopole loop in the renormalized lattice 
unit. For (3 > (3 C the distribution of monopole currents is not uniform as seen in the 
previous subsection. In this case we assume that the distribution is uniform in the 
region restricted to 2R in each spatial direction and Nt in the time direction. Note that 
2R is still larger than Nt near (3 C . Hence the effective renormalized volume which is 
occupied by the longest loop may be given by (2R/n) 3 x (N t /n). l(n) represents the 
range affected by one of the n 3 monopole world line. If the monopole current with a 
range l(n) exists on a link, the nearest current could be put on links more apart than 
l(n). The n 3 extended monopole current with a range l(n) seems to have an exclusion 
volume Z(n) 4 in the renormalized lattice. 

Considering the exclusion volume, the entropy of monopoles may be discussed, in 
the case of n 3 extended monopoles having the range l(n), not on the whole renormalized 
lattice but on a reduced lattice 

(N s /n) 3 x (NJn) , p n n N 

1 ' } m l ' } (for/3 <&), (34) 

where N(l) is an integer depending on l(n). It seems natural to adopt N(l) = [l(n)] + 1, 
where the symbol [l(n)] is the integer not exceeding l{n). Monopoles are seen as if they 
were running only on the reduced lattice. 

l(n) are shown in Figs. 10-12. In the case of l 3 monopole on 16 3 x 4 lattice shown 
in Fig. 10, we see that l(n = 1) is less than 2 for (3 < (3 C and approaches 2 at (3 C . 
l(n = 1) approaches 3 from below as j3 — > j3 c in the case of l 3 monopole on 24 3 x 6 
lattice, whereas l(n = 3) < 1 in the case of 3 3 as seen from Fig. 11. The same regularity 
exists also for any extended monopole on 24 3 x 8 lattice shown in Fig. 12. It is stressed 
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Fig. 10. l(n) versus (3 for l 3 monopoles on 16 3 x 4 lattice. The dotted vertical line denotes f3 c 




Fig. 11. l(n) versus (5 for l 3 (circle) and 3 3 (square) monopoles on 24 3 x 6 lattice. The dotted 
vertical line denotes (3 C . 
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Fig. 12. l(n) versus (3 for l 3 (circle), 2 3 (triangle) and 4 3 (diamond) monopoles on 24 3 x 8 
lattice. The dotted vertical line denotes (3 C . 

that l(n) is always less than or equal to N t /(2n) for j3 < j3 c . This leads us to that 
the reduced lattice always has a time extent 2 = (N t /n)/(N t /(2n)) for (5 < /3 C . Then 
the entropy of a n 3 extended monopole loop for any n should be calculated always 
on a reduced lattice with the time extent 2 for /3 < j3 c . This means that the entropy 
of a n 3 monopole is the same for any n, since the reduced lattice is the same for any 
extendedness n if the original lattice is the same. 

4-6. The entropy of a monopole loop and the deconfinement transition 

In the previous subsection, we have seen that the entropy of a n 3 monopole is the 
same for any n if the original lattice is the same. However, the action is different for 
different extended monopoles at fixed (3. Since the entropy is equal, larger extended 
monopoles are more important to fix the point of the deconfinement transition as seen 
from Fig. 7. 

What is the maximum size of extended monopole? To define dynamical monopole 
currents which are important, the renormalized lattice should have at least two inde- 
pendent lattice sites in the time direction. This means that the maximum size n max 
of dynamical extended monopoles is N t /2. From Figs. 10-12, we see the maximally ex- 
tended monopoles have always l{n max ) < 1 for (3 < f3 c . The reduced lattice is the same 
as the renormalized one for j3 < (3 C in the case of the maximally extended monopole. 

What is the entropy of a long monopole loop on the reduced lattice with N t = 2 
? As shown above, a long monopole loop seems to behave like a random walk with 
a stronger repulsive condition than a simple non-backtracking. However to know the 
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condition correctly is not simple. Moreover, the number of crossing increases as extend- 
edness becomes large. The entropy is very difficult to estimate contrary to the case in 
T = SU(2) QCD. 
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Fig. 13. The histograms of the length of 2 3 monopole loops on 12 3 x 4 lattices. 

Therefore, let us investigate the entropy by a histogram method using Monte 
Carlo simulations. It is easy to get the histogram of the length of monopole loops 
from many configurations generated in the Monte Carlo simulations. The lengths of 2 3 
extended monopole loops are measured on 12 3 x 4 lattice, since they are the maximally 
extended monopole in this case. The histogram is obtained from 3000 configurations 
as shown in Fig.|l^. j3 c is 2.298 on the lattice. The vertical line represent the frequency 
of length. The distribution of the long monopole loops seems to be Gaussian. 

The number of loops of length L through a given point goes asymptotically like 
/i L irrespective of the conditions of random walks.BS Hence the distribution of the long 
loop length may be fitted by the following form 



— aL 2 + bL + c, b = In fi — f\ 



(35) 



where f\ comes from the self energy term in the monopole action and c is a constant. 
The term of —aL 2 should come from all other terms than the self energy and it is very 
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small as seen from Fig. 8. But the term is important because it suppresses proliferation 
of monopole loops in the condensed phase. The linear term may determine the entropy 
ln/x. 

/3 a b 



2.10 1.80 x 10~ 4 2.93 x 10" 1 

2.20 1.41 x 10~ 4 1.76 x 10" 1 

2.30 6.65 x 10" 5 4.93 x 10~ 2 

2.35 5.70 x 10~ 5 2.31 x 10" 2 

Table 3. Least square fitting of the histogram for (3 = 2.10, 2.20, 2.30 and 2.35. 

We plot in Table [3] the values of the parameters of the Gaussian fit to the data 
at (3 = 2.1 ~ 2.35, where the Gaussian shape is seen. We see from the table that the 
linear term at (3 C = 2.298 is already small and it almost vanishes around j3 = 2.35. If 
the balancing of energy and entropy of maximally extended monopole is the origin of 
the deconfinement transition, the predicted critical coupling is (3 C ~ 2.35 which is near 
to the actual critical coupling (3 C = 2.298. Considering the finite size effects near /3 C , 
the agreement is fairly good. 

From the fit, one can estimate roughly the value of the entropy using the values 
of /i(/3), i.e., A (2.35) ~ 1.3 and f x (2.298) ~ 1.16. We get 1.16 < In// < 1.3. 

In Fig. 14, we plot f\ for maximally extended monopoles on various lattices in 
comparison with the entropy estimated above. We find that the balancing of the en- 
ergy and the entropy of the monopoles can explain the transition points of all lattices 
approximately. There are, however, many things to be studied. Nevertheless, it is very 
interesting that the deconfinement transition can be understood roughly by the simple 
idea of balancing of energy and entropy of monopole loops in T ^ SU(2) QCD. 

§5. Conclusion and remarks 

1. In T y^ QCD, it is found that a long monopole loop exists in each configuration 
and all other loops are short in the confinement phase. They distribute almost 
uniformly in the confinement phase. Table 1 shows that the number of short loops 
decrease in the case of 2 3 monopoles in comparison with the case of l 3 monopoles. 
To study the continuum limit, we have to make the time extent N t large. Then 
large extended monopoles are important and short loops would disappear rapidly. 
It is expected that long monopole loops alone are responsible for the confinement 
mechanism of quarks. This is consistent with the result that the string tension is 
reproduced only by long monopole loops.cPW 

2. Almost all magnetic charge is ±1 and the loop length is rather restricted in 
comparison with all link numbers of lattice. There appears to work some repulsive 
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Fig. 14. Coupling constant /1 for various extended monopoles. The dotted vertical lines denote 
(3 C for Nt = 4, 6 and 8 from the left. 

force between currents, so that monopoles may have an exclusion volume around 
them. Considering the exclusion volume determined from the distribution data, 
we have derived that the entropy of extended monopoles may be independent of 
the extendedness if the original lattice is the same. 

3. InT^ QCD there exists a maximally extended ((A^/2) 3 ) monopole contrary 
to the case of T = SU(2) QCD where larger and larger extended monopoles 
become important. The critical temperature becomes necessarily finite. Balancing 
of energy and entropy of the maximally extended monopoles can explain roughly 
the deconfinement transition of various lattices. However more intensive studies 
are needed to clarify the transition mechanism definitely. Especially, the finite 
size effect near f3 c should be clarified in details. 
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